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Abstract. In this paper we study the asymptotic behaviour of solutions of a system 
of N partial differential equations. When N = 1 the equation reduces to the Burgers 
equation and was studied by Hopf. We consider both the inviscid and viscous case and 
show a new feature in the asymptotic behaviour. 
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1. Introduction 

We consider the system of generalized Burgers equations for N unknown variables u = 
• • -,Un), 

v 

( u j)t *r(c, U)(Uj) x = ^ {uj)xxi j = 1 ;2,... ,N, (1-1) 

where c = (c\,C 2 , ■ ■ ■ ,cn) is a constant vector in R N and cj(c,m) = c kUk is the usual 
inner product in R N . We study the solution of (1.1) with initial conditions 

uj(x,0) =u 0j {x), j = 1,2,... (1.2) 

When v>0, (1.1) is a system of nonlinear parabolic equation describing the interplay 
between nonlinearity and diffusion, v being the viscosity parameter. 

When v = 0, the system (1.1) is hyperbolic with coinciding wave speeds ct(c,m), and 
the nonlinearity and the nonconservative form makes the initial value problem complex. 
Indeed the system being nonlinear, solution cannot be continued as a smooth solution even 
when the initial data is smooth. Further if N > 1, the product cr(c, u) • ( Uj) x is nonconser¬ 
vative and does not make sense in the usual distributional sense. The solution should be 
understood in a generalized sense. 

For the case v > 0, Joseph (5) used a generalized Hopf-Cole transformation to linearize 
the system of equations (1.1) and solve (explicitly) with the conditions (1.2) in terms of a 
family of probability measures du^ (v). These measures depend on the initial data (1.2) 
in a nonlinear and nonlocal manner and takes the form 








(1.3) 
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where 


I(y)= [ cr(c,uo(z))dz. 
Jo 


(1.4) 


When v > 0, the solution of (1.1) and (1.2) was shown to be 



(1.5) 


Let nj(x,f) be the solution of (1.1) and (1.2) given by (1.5). It was shown in 0 that 
when UQj is Lipschitz continuous, for each t > 0, except for a countable x the limits 


u;(x,t) = lim u v Ax,t ) 
v-*o J 

exist and is given by the formula 


uj(x,t) =u 0 j(y(x,t)), j=l,2,...,N 


( 1 . 6 ) 


where y(x,t ) is a minimizer of 


min (/(y) + (x — y) 2 /2t) 

oo<V<oo 


(1.7) 


and I(x) is given by (1.4). 

Using some ideas from the earlier works of Joseph Qj) and LeFloch Uj(x,t)j=\. 2 .....N 
was shown to be a solution of an inviscid case (v = 0) in (1.1) with initial data from (1.2), 
the nonconservative product was justified in the sense of Volpert product £02 (see Dal 
Maso, LeFloch and Murat 0 for a generalization of Volpert product). In 0 solution for 
general initial data is constructed in the sense of Colombeau (T|. 

The aim of the present note is to study the asymptotic behaviour of the solution for 
the parabolic (viscous) case as well as the hyperbolic (inviscid) case. Study of asymptotic 
behaviour of solutions is important, see 0 and the references therein for the parabolic 
case and jTJ for the inviscid case. When N = 1, this system is the celebrated Burgers equa¬ 
tion and explicit solution and its asymptotic behavior as t tends to infinity and diffusion 
parameter v —> 0 was studied by Hopf |3j. We show that Hopf’s analysis give asymptotic 
form of the solution for the viscous case. When N = 1 and c i / 0, it is well-known from 
the work of Lax (7) that for solution of inviscid Burgers equation with initial data sup¬ 
ported in the compact interval [—£,£], £ > 0 the solution decays at the rate 0(t~?) and 
support spreads at a rate 0(t ?) for large time. From an explicit computation we will show 
that the decay rate is not true in general for the present case, but still the support spread at 
the same rate. We start with the viscous case. 


2. Asymptotic behaviour with viscous term 

In this section we study the asymptotic behavior of solution of (1.1) and (1.2) when v > 0 
and fixed. On the initial conditions UQj{x) 1 j = 1,2,... ,N, assume that =/(°°), 

lim^_<>o/(x) =/(—°°), Xim x ^ooUQj(x) = uoj(°°), lim^_ooMo ; (x) = uqj (—°°) exists and is 
finite. With these assumptions, we shall prove the following result. 


Asymptotic behaviour of solutions 


511 


Theorem 2.1. The solution u v (x,t ) = {u\(x,t),U 2 (x,t),... ,u^(x,t)) of ( 1.1) and (1.2) 
has the following asymptotic behaviour as t tends to infinity : 


u V :(x,t ): 


/ \ r x/\T(tv) _ 
MO; He v i_oc e 


y 2 j 

Tdv- 


MO 


oo e 




s/y/ltv) ‘ 


_v 2 

Tdv 


-2(°°) r x l\fi tv ) y2 j 

e~ —f_J e t dy - 


-n—) 


roo 

J v/vH) e 


V 2 , 

T dv 


( 2 . 1 ) 


Proof First we note that the solution u?(x,t),j = 1,2,... ,1V of (1.1) and (1.2) is given by 
(1.5) where the measure djU^. ( ,(y) is given by eqs (1.3) and (1.4). Writing explicitly the 
formula, we get 




f-ooU 0 j(y)e 


_ v /(y)+ 


(*-y) 2 


I dy 






J dy 


Setting <; = x/s/vt, and then making a change of variable z 
as y, we get 


\Jvt 


and renaming z 


'(x,t) = 


f-coUoj(Vvt(Z-y)e 


\ i(Vn(S-y } +y 2 


7 2 


dy 


r /(Vvt(j-y)) | 2 / 2 l 

n.e L V +>/ 2 J d y 


( 2 . 2 ) 


Now we split the integral in (2.2) in the following manner: 


r 00 __ r /(Vv7(g-v)) 2 /o] 

/ M 0 ;(\/v7(Hy)) e ^ v ' ^dy 

J —oo 


HVvHi-y)) | ..2 /2 

-V))e l - +y J dy 


I u 0j {Vvi{^-y)) 

J —oo 

ooo 

+ U 0 j(Wt(!;-y)) 

Jq+o 

r£+8 , r /(y^7(4-y)) 2 /2 ] 

+ / M 0 ;(Vv7(Hy))e ^ v ‘ J dy. 

J £—8 


7(\/v?(g-y)) | 2 / V 

— v))e L v +y/ Jdy 


(2.3) 


Now we fix 5 > 0 and study each of these integrals as t tends to infinity. We have under 
the assumptions of the theorem on I(x) and uq j{x), as t tends to infinity: 


u 0j (Vvi{Z-y)) 
f « 0j «-y))e 

J £,+8 


_ r nvvta-y)) , v 2 /2 i /(+«,) r 

l v +y/2 \d [ywe“ v k 0j .(oo)/ 

r ; (^-y)) 2 /2 l 

L v J dy«e v U(\A — 


«"* .,2/ 


lim sup 


/■€+S _ r /(y^7(^-y)) | 2/ 2 ] 

/ m 0 ;(v / w(^ - y))e ^ v "Jdy 

J £—8 


e -y ( 2 dy. 
MO;(-°°) [ e“ r/2 dy, 

= 0(5)- 
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Now let t tend to infinity and then 8 tend to 0 in (2.3). We get 

u 0j {Vw{Z -yVe-l^tt-yVv+yV 2 ]^ ~ e -^u 0 j{°°) J" e'^dy 


+ e £ \q J e y / 2 dy. 


(2.4) 


Similarly 

£ 


e _[5ViM +)?/2 ] dv!se=M fi 


J e ^'dv 4- e i 1 J e ^ 2 dy. (2.5) 


We observe that due to our assumption on 7(x), this limit in (2.5) is a postive real number 
and hence letting 1 tend to infinity in (2.2) and using (2.4) and (2.5) we get the result (2.1). 
The proof of the theorem is complete. □ 


Remark. An interesting case here is when the initial data uqjJ = 1,2,...,A satisfies 
the following conditions. no ; (°°), «o/(—°°) are nonzero and there is a cancellation in 
<j(c,uq)(x) = Yii c k ll ok(x) so that this quantity is integrable. 


3. Asymptotic behaviour of solutions of generalized Hopf equation 

In this section we study the asymptotic behaviour of vanishing viscosity solutions of 

(uj), + o(c 1 u)(u j ) x = 0 ( 3 . 1 ) 

with initial data 

Uj(x,0) = uoj(x), ( 3 . 2 ) 

for j = 1 , 2 , ...,1V. We recall the definition of cr(c,n) namely <j(c,u) = c A M /t where 
c = (ci,C2,... ,cn ) a given constant vector and n = (u\,U 2 , ■ ■ ■ ,un )■ the unknown vector 
variable. 

It is easy to see that u y (x,t),j = 1 , 2 ,... ,Ais a solution of ( 1 . 1 ) and ( 1 . 2 ) iff <7 v (x,f) : = 
<j(c,u v (x,t)) satisfies the Burgers equation 



with initial condition 

cr(x,0) = cr(c,no)(x). 

By Hopf-Cole transformation, the solution can be written in the form 

<y v {x,t) = J <7(c,Mo(y))d^ ) (y)- 

As in j6j it is easy to see that, when uqj is Lipschitz continuous, for each t > 0, except for 
a countable number of points of x, the limit 

cr(x,f) = lim o v (x,t ) 
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exists and is given by the formula 

o{x,t) = o(c,u 0 (y{x,t))) = c 
where y(x,t ) is a minimizer of 

min (I(y) + {x — y) 2 /2t) 

—oo<y<oo 

and I(x) is given by (1.4). Here we remark that this formula is slightly different from that 
of Hopf 0. Once a{c,u) is known, (3.1) can be treated as N scalar linear equation with 
discontinuous coefficient. For any function a(x,t) on [(x,f): —°°<x<°°,t >0], ||a(.,f)||<«, 
denotes the essential supremum of a(x,t) with respect to the space variable x keeping the 
time variable t > 0 fixed. We define the left boundary curve x = s' (/) and right boundary 
curve x =s + (t) of the support of a{x,t) where s~(t) = sup[y: a(x,t ) = 0 for all x < y] and 
s + (t) = inf[y: a(x,t) =0 for allx > y]. It is well-known from (7]| that when the initial data 
is supported in [—/,/], being solution of the inviscid Burgers equation, (7(c. u) have the 
following estimates. There exist constants A > 0 and C > 0 which depend on / and || mo ||°o, 
such that 

- 1-At? <s~{t) <s + {t) <l+At?, 

IK.,/)||oo <CfT. (3.3) 

We shall prove that for u/,j = 1,2,... ,N again the same estimate holds for the spread of 
support but the decay result is not valid in general. First we have the following theorem. 

Theorem 3.1. Let Uj{x,t),j = 1,2,... ,N be the vanishing viscosity solution of (3.1) and 
(3.2) with initial data supported in [—/,/]. Let x = sj (f) and x = sj (t) are the support 
curves for uj(x,t). Then there exists a constant A > 0 which depend on l and ||mo||oo, so 
that the following estimate holds for t 1 , 

-l —At? < sj (/) < sj(t) < l +At 2 . 

Proof. The proof easily follows from the fact that the characteristic speed < 7 ( 7 :, u) = 0 
outside the region —l—At? < x < l +At T So in this region the characteristics connecting 
the points (xf) to a base point (y, 0 ) are parallel to the f-axis and hence is of the form 
x = y and since the solution is constant along the characteristics, we have u(xf)=uo{y) = 
U()(x). When {xf) is outside —l—At? < x < l +At?, the x co-ordinate of the base point 
of the characteristic, y, lies outside [—/,/] where uo(y') is zero and hence u(x,t) is zero. 
The proof of the theorem is complete. □ 

Next we show that solution does not decay in general, by giving an example. Here we 
construct vanishing viscosity solution of (3.1) with initial data of the special form which 
is supported in a compact set, namely 

{ 0 , if x<—l 

UQj , if — / < x < l (3.4) 

0 , ifx>/ 

where UQj is a constant and / is a positive real number. Let wj(x,f). 7 = 1 , 2 ,...,at be the solution 
of (1.1) with the initial data (3.4). We shall prove the following. 
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following form: 
When Oo = 0, 


When Oq < 0, 


uj(xf) = < 


- 2-*= 

_\CjUQj, then Uj\ 

0, 

if x < — / 

M0;, 

if — / < x < / 

0, 

if x > / 

0, 

if x < ?j- • 

M0;, 

if - / 

«;7 
°0 ' 

if Ool < x 

o, 

if x < / — 

“0; 

°o 

if/-(-4 

.0, 

if x > / 


(3.5) 


- 4 / 

Co 


-41 

°o 


-41 

°o 


- 4 / 

CTO 


(3.6) 


-41 

Co 


When Oq > 0, 


lim u v Ax,t) = ^ 


v—>o 


if x > ^ 


if o 0 r - i 


41 


°o 

T‘ 1 l ’ f < 5b 

^.£±i, if -l <x<a v-l,t<M. 

0, ifx>/ + (4/o 0 /)V > ^ 

if —l<x<l + (4/ooO^f > 
if x < —/ 


(3.7) 


v0, 


Proof. To prove (3.5)-(3.7) we use the formula (1.5) to get explicit solution of (1.1) and 
(3.4) in the form 


u] 0,0 = 


UjO [^2f L+o ' 0 - v ] dy 


5-L e-^dy + /!, e“ v t^+w] d y + j; 

which can be written in the form 

UjO f-/ e ~^ [^f L+CT ° ) '] dy 


(x-y ) 2 , 

/ e ~^sv~ d y 


(3.8) 


uj(x,t) = 


1 r i 1 f h-/) 2 I q .,,.1 , ./ h-h 2 

z+f ,e v L 3? dv -I- P ,e ' 


(3.9) 


(2f7rv)z + /_,e 


ldy + /_ ; e "Tfv - dy 


To study the limit, we rewrite this formula in a convenient form by introducing the func¬ 
tions 


A l,o 0 ( x A) = 


, .1 °n < °o* / fOo —X—/ \ 

,f) = (2 tv) 1 e^y v erfc I -—— ] 


(3.10) 
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and 


where 


V , s / . l a n‘ apt , /fCo —X + /\ 

B l<j 0 ( x f) = (2fv)?e^ v er / c - , 

V ( 2 #v )2 y 


er/c(y) = [ e '"dy. 
Jy 

We can rewrite (3.9) as 

t/y (x,f) = 


u Oj( A lc 0 ( x ’t) ~ 


{2t7tv)i +Aj (Jo (x,t)-Bj OQ (x,t) +Aj 0 (x,t)-Bj 0 (x,t) 
Using the asymptotic expansions of the erfc , namely. 


er/c(y) = 


and 


1 1 

2y ~ 4y3 


1 


1 


e r/c(-y) = W^^--^3+o^ 
in (3.10) and (3.11) we have the following as v —> 0: 


(tv) 


X+Oq/ 


e 2vF 


A la 0 { x f) '■ 


i °o f _ £o£ 


(WV) 2e ^r 


if — / — x + (Tot > 0 
if — l — x + Got = 0 • 


2 

if -l-x+o Q t < 0 


M 

/—X+(jQ t 


e 2 v? 


i <V _ 2a£ 


B l<j 0 ( x f) ~ < (^)le^- 


if / - x + Cof > 0 
if / — x + Oo f = 0 • 


I On* oirtjc v -2 

(2 7 Tfv)^e^v~~ v + /_ v+ g f e~ , if/—x+(7o f <0 


(3.11) 


(3.12) 


(3.13) 


(3.14) 


(3.15) 


-jr+oo? 

It is straightforward to check the formulas (3.5)-(3.7) using (3.14) and (3.15) in (3.13). 
When cjq = 0, (3.13) becomes 


uVxf) = 


u Oj( A lo( x f) - B lo( x f)) 

(2 tnv)? 


(3.16) 


Now take the region x < —/. Then — l — x > 0 and / — x > 0, and using (3.14) and (3.15) 
in (3.16) we get 


2 2 
tv iL tV 

V/ ^, —/—x l—X 

Uj [x, t ) « UQj - 7 - 


(2 otv)2 


and hence we have 


lim u v Axf) = 0, x < —/. 
—.0 •* 
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In the region —l < x < l, — / — x < 0,1 — x > 0, 


(: 2mv )2 


tV ■ fv ■ 

—V— Q^Itv — 

—l—x l—x 


u V j(x,t ) w UQj 

J (2ntv)i 

Thus we get 

lim u v Ax,t) = UQj, —l<x<l. 
v->o J J 

In the region x > L we have — / — x <0,1 — x <0 and 


;(x,t) ~ UQj 


(2 ntvY- + - (2mv)^ - ^ew 

(2 7ttv)y 


and hence it follows that 


lim u V :(x,t) = 0, x > l. 
v—>o J 

This completes the proof of the theorem for Oo = 0. The case <7 / 0 is similar and is 
omitted. □ 


Remark. Thus if we take the initial data (3.4) which has compact support the decay of 
the vanishing viscosity solution uj(x,t),j = 1,2,... ,N depends on the initial speed Co = 
Yl\ c^UQk- If tTo = 0, the vanishing viscosity solution does not decay. Indeed 

sup \uj(x,t)\ = \u 0 j\. 
xeR 1 

On the other hand, from the above theorem it follows that for the case (To £ 0, the solution 
decays, namely 

sup \uj(x,t)\ = o(t~~ J. 
xeR' V y 

We conclude with a remark on the solution for (3.1) with the Riemann type initial data 


(■ Uj)(x, 0 ) 


ujl, ifx<0 
ujr, ifx>0, 


(3.17) 


where uji and ujr are constants for / = 1 . 2 .... ,n. 

Let (T/ = Ylk= \ c k u ji. and Or = D" i c kUkR- Then the vanishing viscosity solution for 
the Riemann problem (3.1) and (3.17) takes the following form 0- 
When Ol< Or, 


UjL, 

if x < Ojj 

UjR-UjL X I u jL^R-Uj R a L 

a R -a L t o R -a L ’ 

if opt < x < ORt 

UjR, 

if x > ORt 
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When ol = Or = a. 


Uj(pC,t) 

When ai > a r, 


uj(x,t) 


u jL, 

if x < a t 

UjR, 

ifv > a-t. 

UjL, 

if x < ° L +° R -t 

UjL+UjR .f ^ _ o l +Cr t 

UjR , 

ifx> ° L + CR -t 
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